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Abstract 



In formal scattering theory, Green functions are obtained as solutions of a dis- 

^ ■ tributional equation. In this paper, we use the Sturm-Liouville theory to compute 

\G • Green functions within a rigorous mathematical theory. We shall show that both the 

Sturm-Liouville theory and the formal treatment yield the same Green functions. We 

f~^ ' shall also show how the analyticity of the Green functions as functions of the energy 

^D ■ keeps track of the so-called "incoming" and "outgoing" boundary conditions. 

PACS number(s): 03.65.Nk, 23.60. +e* 

^ '■ 1 Introduction 

Green functions are essential tools in mathematical physics |l|, . Most books on scattering 
theory ^, Q only provide a computational procedure to obtain Green functions without any 
reference to a mathematical setting. Here we provide that mathematical framework for the 
example of a square barrier potential. In this example, the mathematical framework is given 
F^. ! by the Sturm-Liouville theory 0. 

We consider a square barrier potential of height Vq, 






<r < a 


Vo 


a <r <b 





b < r < oo 



V{r) = < Vo a<r <b (1.1) 

[ 6 < r < oo. 

Due to the spherical symmetry of the potential, spherical coordinates come in handy. The 
expression for the Hamiltonian in spherical coordinates (for s-waves) is 

h^-f-^+V{r). (1.2) 
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In order to obtain a linear operator from the formal differential operator ( |1.2| ), we need to 
define a domain in the Hilbert space -^^([0, oo), dr) on which h acts. The domain we choose 
is 

V{H) = {fir) I {hf)ir), fir) G L'i[0, oo), dr), /(O) = 0, /(r) G ^^^[0, oo)} . (1.3) 



On this domain, the formal differential operator ( p.. 21) induces a self-adjoint operator H 
(cf. [|]). The spectrum of this operator is the positive real line [0, oo) (cf. 0). 

Both formal and Sturm-Liouville approaches use the solutions of the time-independent 
Schrodinger equation, 

dr'^ 

as their basic source of information. The difference is that physicists, fearful of using complex 
energies when working with self-adjoint Hamiltonians, solve ( |1.4|) using real energies. The 
Sturm-Liouville approach uses complex energies. This is in no contradiction with the self- 
adjointness of the Hamiltonian, since the eigenfunctions that correspond to complex energies 



-— + V{r))a{r,E) = Ea{r,E), (1.4) 



lie outside the domain (|1.3| ) on which the Hamiltonian is self-adjoint. 

In the next section, we show that the formal and the Sturm-Liouville approaches yield 
the same Green functions. In our calculations, we will use the following branch of the square 
root function: 

y^ -.{E eC\ -n < arg(E) < vr} i — >{E eC\ -n/2< aig{E) < n/2} . (1.5) 

2 Computation of the Green Functions 

2.1 Formal Approach 

To compute the (radial) Green functions, physicists solve the following distributional equa- 
tion: 

d' 



^^2 + ^(^) - ^ ) G{r, s- E) = -6{r - s) , E G [0, oo) , (2.1) 

subject to certain boundary conditions. Equation (^TT| ) says that for r ^ s, G{r, s; E) obeys 
the radial Schrodinger equation ( |1.4D and the following boundary conditions: 

G{a-,s;E) = G{a+,s;E) (2.2) 

G{b-,s;E) = G{b+,s;E) (2.3) 

-^Gia-,s;E) = ^G{a+,s-E) (2.4) 

or or 

^G{b-,s;E) = ^G{b+,s;E). (2.5) 

or or 

At r = s it is continuous, but its derivative has a discontinuity of 1, 

—G{s + 0, s; E) - —Gis - 0, s; E) = 1 , (2.6) 

or or 



due to the delta function. At the origin, the Green function must be regular, 

G{0,s;E) = 0. 



(2.7) 



There are two linearly independent Green functions that we are mostly interested in: 
the incoming and the outgoing Green functions. The outgoing Green function G^ satisfies 
Eq. ( p.l|) , the boundary conditions ( |2.2| )-( p77D , and an "outgoing boundary condition" at 
infinity, 

G+{r,s;E)^e'^\ for r -^ oo . (2.8) 



The boundary condition ( p.8| ) means that very far from the potential region the function 
G^ behaves as an outgoing wave. The expression for this outgoing Green function is f^ 



G^{r,s;E) 



x(r<;E)^+(r>;E) 



(2.9) 



where r< and r> refer to the smaller and to the bigger of r and s, respectively, x{f"i E) is 
the solution of the Schrodinger equation that vanishes at the origin, f2+(r; E) is the solution 
of the Schrodinger equation that satisfies an "outgoing boundary condition" at infinity, and 
W{Xy ^+) is the the Wronskian of x and ^2+. More precisely, x{^] E) satisfies Eq. ( |1.4|) and 
the boundary conditions 

(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 

Its expression is given by 
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X(r; E) = { Ji(E)eV(^-^o)r ^ J^^^•^^^^^iE-Vo)r 
J^{E)e'^'- + ME)e-'^'- 



< r < a 

a < r < b 
b < r < oo . 



(2.15) 



The functions J\-J^ are such that x(^; -^) satisfies the boundary conditions ( p.lO| )-( ^yi^ ) 
and are given in the Appendix. The function f2+(r; E^ satisfies Eq. ( |1.4| ) and the boundary 
conditions 



fi+(a+;^) 

fiV(a+;^) 

fi+(6+ : E) 

n'^ib+;E) 

n+{r;E) 



Q+{a-;E) 
n'^{a-;E) 

n+{b--E) 

n'^{b- : E) 



OO 



(2.16) 
(2.17) 
(2.18) 
(2.19) 

(2.20) 



The solution fi+(r; E) reads (see the Appendix for the exphcit expressions of Af-Al) 

( Ai{E)e' ^' + A t{E)e-'^' < r < a 

fi+(r; E) = i ^+(E)eV(^-^o)r ^ ^+(^)e-V(^-^o)r ^ < ^ < 5 (2.21) 

[ e'^'' b<r < 00. 

The Wronskian of x(^; -S) and f2+(r; £") is 

W{x, fi+) = 2ty/EME) . (2.22) 

The outgoing Green function is therefore given by 



The incoming Green function G satisfies Eq. (|2.1| ), the boundary conditions (p.2|) -i 
and an "incoming boundary condition" at infinity, 

G- (r, s; E) ~ e"^^" , for r ^ 00 . (2.24) 

Condition ( |2.24| ) says that far away from the potential region G^ behaves as an incoming 
wave. The expression for this incoming Green function is 

G (r, ., E) ^(^^^_) , (2.25) 

where x('"j E) is the solution of the Schrodinger equation that vanishes at the origin, 
VL_{r\E) is the solution of the Schrodinger equation that satisfies an "incoming boundary 
condition" at infinity, and W{Xi ^-) is the the Wronskian of x a-^id VL_. The eigenfunction 
x{r] E) is given by (|2.15|) . The eigenfunction fl_{r]E) satisfies the Schrodinger equation 



subject to the following boundary conditions: 

n^{a+;E) = fi_(a-;E) (2.26) 

n'_{a+;E) = n'_{a-;E) (2.27) 

n-{b+:E) = n^{b-;E) (2.28) 

^'_{b+;E) = ni{b- : E) (2.29) 

Q_{r;E) ~ e''^' , r ^ 00 , (2.30) 

and its expression is given by 

A^ {E)e'^'' + A2 (E)e-^^^ < r < a 

Q.{r;E) = <J ^3 (E)eV(^-^o)r ^ ^j(^)e-V(^-^o)r a<r<b (2.31) 

a-iy/Er b < r < 00 . 



The functions Ai-A^^ make Q^{r; E) satisfy the boundary conditions ( p.26| )-( p.30|) . The 
Wronskian of x(^; E) and r2_(r; E) is 

W{x, n_) = -2iy/EJ^{E) . (2.32) 

Therefore, the expression of the incoming Green function is 

2.2 Sturm-Liouville Approach 

In this section, we compute the Green functions using the Sturm-Liouville theory. The Green 
function appears as a kernel of integration when we write the resolvent of the Hamiltonian 
H as an integral operator, 

/•oo 

{E-Hy^f{r)= G{r,s-E)f{s)ds, E ^ [0, oo) . (2.34) 

Jo 



From Eq. ( |2.34| ) one can formally derive Eq. ( |2.1| ). 



As mentioned in the Introduction, physicists avoid computing Green functions using 
complex energies, but rather impose the incoming and outgoing boundary conditions. The 
Sturm-Liouville theory deals directly with complex energies in a way that keep track of the 
"in" and "out" boundary conditions. 

The procedure to compute the Green function for our Hamiltonian H is given by the 
following theorem (see Theorem XIII. 3. 16 in [^): 

Theorem 1 Let H be the self-adjoint operator derived from the real formal differential 
operator (|1.2|) by the imposition of the boundary condition /(O) = 0. Let Im(ii^) ^ 0. Then 
there is exactly one solution x(r; E) of {h — E)a = square-integrable at and satisfying 
the boundary condition /(O) = 0, and exactly one solution Q{r] E) of {h — E)a = square- 
integrable at infinity. The resolvent {E — H)~^ is an integral operator whose kernel (j(r, s\ E) 
is given by 

{ x{r;E)nis;E) 
w(x,n) 
(2.35) 
x{s;E)n{r;E) ^ ^ 

w{x,n) ' ^ b- 

The incoming and outgoing Green functions are defined by 

G^{r,s;E) = \imG{r,s]E±iii). (2.36) 



From Eq. ( p.36| ) it follows that the Green functions G^ are just the boundary values on 
the real axis of the kernel (|2.35|) , which is a function of complex variable. Thus G^ keep 



track of what is "incoming" and of what is "outgoing" by specifying which side of the cut 



of the resolvent we are on. In our example, the cut of the resolvent, i.e., the spectrum of 
the Hamiltonian, is the positive real axis. 

First we compute the kernel (|2.35|) in the region Im(i^^) > 0. In this region, the eigen- 
function x(r; E) in (|2.35| ) satisfies the boundary conditions 



X{0-E) = 0, (2.37) 

xir;E)EAC\[0,oo)), (2.38) 

x{r; E) is square integrable at . (2.39) 



The boundary condition (|2.38|) implies the boundary conditions (|2.11|) - (|2.14|) , and ( p.39| ) is 



automatically fulfilled. This means that the eigenfunction xi.f'^E) of Theorem 1, which is 
unique, is the same as the eigenfunction of Eq. ( |2.15|) , although now E is allowed to run 
over the upper half-plane of complex numbers. 

The eigenfunction VL{r\E) of Theorem 1, that we denote by VLj^{r]E) for the case of 
Im(£') > 0, satisfies the Schrodinger equation and the boundary conditions 

n+{r- E)eAC\[Q, oo)), (2.40) 

fi+(r; E) is square integrable at cxd . (2-41) 



The boundary condition ( p. 401) implies the continuity conditions ( p.l6|) - (p.l9|) . Condition 



(|2.41| ) is, for Im(i?) > 0, equivalent to ( p.20| ). Thus the (unique) function VL^{r]E) of 
Theorem 1 coincides with (|2.21|) , E being now any complex number with positive imaginary 
part. Thus, for Im(i5) > 0, the expression of G(r, s; E) in ( |2.35| ) is given by ( |2.9| ), although 
now E is allowed to be any complex number in the upper half-plane. 

For \m.{E) < 0, the situation is similar. The eigenfunction x(^; E) of Theorem 1 is also 
given by (|2.15|) . The other eigenfunction in ( p.35| ), that we denote by VL-{r\E), coincides 



with the eigenfunction of Eq. ( |2.31[ ). Needless to say, E can now be any complex number 



in the lower half-plane. Therefore, the kernel (|2.35|) is given by ( |2.33| ), E now being any 
complex number with negative imaginary part. 



If we now compute the limits in ( p.36|) , we obtain exactly the same outgoing and incoming 



Green functions as those obtained in the previous section. 

3 Conclusions 

We have seen that the Green functions obtained by applying the Sturm-Liouville theory 
are the same as those used in formal scattering theory. We have also seen that in order to 
find out what is "incoming" and what is "outgoing" , the Sturm-Liouville theory substitutes 
the standard physicists boundary conditions by statements on the analyticity of the Green 
functions with respect to the energy variable — knowing which side of the cut we are on is 
tantamount to knowing what is "incoming" and what is "outgoing" . 

These conclusions are not restricted to the square barrier potential. Actually, they hold 
for potentials that decrease fast enough at the origin and at infinity, and that do not have 
too many discontinuities. 

6 



Acknowledgments 

The author wishes to express his gratitude to M. Mithaiwala for proofreading the paper and 
making many invaluable suggestions, and to Prof. J. P. Antoine and Dr. T. Kuna for intro- 
ducing him to the Sturm-Liouville theory. Financial support from the Welch Foundation is 
gratefully acknowledged. 



Appendix 



The functions that appear in the expressions of the eigenf unctions in Section |2.1| are given 
by 



J,{E) = ^e-'^" 



Ve \ 

— = cos(vEa) I 

^cos(VEa) 

1 - VW^\ ^-^^fi^^)^j^^E) 



ME) = L^Vi^^)'^ (siniVEa) - -, ^_^ 
2 \ i^/{E-Vo 



ME) = ^e 



iVEb 



^ _ VWZM] e^v^ 



Vo 



E 



^'Ji{E) 



l+v{(^M)e-xA^^''j, 



{E) 



(3.1) 
(3.2) 



(3.3) 



(3.4) 



At{E) 

At{E) 
AUE) 



2 I ^{E^V^) 



E 



2 I ^{E - V,) 



i^/Eb 



iVEb 



-e 
2 



i^/Ea 



^^ V(E-K))\ ,y(^Z^ 



(E) 



^ _ V{E-Vo) \ ^_,^/lE^^)aj^+ 



(E) 



At{E) 



-e 
2 



:VEa 



\fE 



\E) 



(3.5) 
(3.6) 



(3.7) 



1 + ^^^^^) e~^Vi^^)^At{E) 



(3. 



Al{E) 
A-,{E) 



2 I ^{E - V, 



2 I y/{E - Vo 



-i^/Eb 



-iy/Eb 



-iy^a 



A^{E) 



AvEa 



1 + v{(^M j e^Vm^)'^A^(E) 
h (l - ^^^^^) e-^^/^^^^^A-AE) 
1 - V^^M^j e^V(^^)-AAE) 
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